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The nonlinear algo rithms proposed re cently by Abrams 
and Lloyd [Report No. quant-ph/9801041] are fast but make 
an explicit use of an arbitrarily fast unphysical transfer of in- 
formation within a quantum computer. It is shown that there 
exists a simplification of the second Abrams-Lloyd algorithm 
which eliminates the unphysical effect but keeps the algorithm 
fast. 



I. INTRODUCTION 

Any systematic procedure associating with an arbi- 
trary number iq . . An-i, = or 1, another number 
/(?o ■ • -^n-i) can be termed an algorithm. The idea of 
quantum computation rests on the observation that a bi- 
nary number zo . . . in-i can be represented by a vector (a 
qubinary number) . . . \in~i) representing an uncorre- 
lated state of n distinguishable two-level quantum sys- 
tems. A quantum algorithm is essentially an algorithm 
based on a qubinary representation of numbers j^J^ . 

Typical quantum algorithms use unitary operations 
and projections . This is motivated by the unitarity 
of the standard Schrodinger dynamics and the so-called 
projection postulate. The latter postulate is typical of 
the Copenhagen interpretation of quantum mechanics 
and is not essential to quantum computation (quantum 
algorithms would look different in, say, the Many Worlds 
interpretation of quantum mechanics but their funda- 
mental properties would not change). The Schrodinger 
(linear and unitary) dynamics is not the only dynamics 
one encounters in quantum theories. Dynamics in the 
Heisenberg picture is typically nonlinear (in the sense 
that operators depend nonlinearly on initial conditions). 
Typical effective dynamics of quantum optical systems 
(such as two-level atoms) is irreversible and nonunitary. 
There are also many situations where states of quantum 
systems evolve in an effectively nonlinear way (optical 
solitons, Hartree-type approximations). Finally, there is 
still no proof that the present-day quantum mechanics 
is not an approximation to a more exact nonlinear the- 
ory, and various versions of such a nonlinear generaliza- 
tion have been proposed. It should be stressed that the 
popular oppinion stating that all nonlinear extensions of 
quantum mechanics must lead to logical absurdities does 
not find support in a detailed analysis of nonlinear "no- 
go" theorems (for a brief discussion cf. j|]). 

The assumption of fundamental quantum linearity can- 
not be regarded as a consquence of experimental data be- 
cause it is quite typical that a consistent theory is prior to 
experiments. This point was clear to Wigner H who was 



simultaneously one of the first to associate fundamen- 
tally nonlinear phenomena with a theory of brain func- 
tioning [Q. Since there is no doubt that human brain is 
a physical system, there is almost no doubt that at least 
some of its aspects have to be described by quantum 
mechanics. The idea of qubinary mathematics (includ- 
ing diferentiation and integration) can be traced back to 
Orlov's works on a "wave logic" of conciousness Q . On 
the other hand, there are serious arguments of Penrose 
|p|,p^ for non-algorithmic (in the classical sense) ingre- 
dients in brain activity. Human and animal brains are 
systems that seem to possess a feedback-type property 
of self-observation, and feedback effects are typically as- 
sociated with a nonlinear evolution. 

Although the above problems may appear somewhat 
far from standard quantum physics, they naturally lead 
to the question of possible consequences of a nonlinear 
quantum dynamics for the theory of quantum compu- 
tation. Theproblem was recently addressed by Abrams 
and Lloyd |ll|| who showed that a nonlinear evolution in a 
Hilbert space of states of a quantum computer leads nat- 
urally to polynomial-time solutions of NP and #P prob- 
lems. The Abrams-Lloyd argument was based on a gen- 
eral property of nonlinear evolutions in Hilbert spaces, 
namely the non-conservation of scalar products between 
nonlinearly evolving solutions of a nonlinear Schrodinger 
equation. This effect (in the literature called a mobility 
phenomenon) was discussed in great detail by Mielnik 
in his analysis of nonlinear motion semigroups (com- 
pare also ]l3t). 

The mobility effect was used in in two algorithms. 
In the first of them the Authors chose a nonlinear evolu- 
tion which has |0) as a fixed point, but any superposition 
of |0) and |1) transports towards |1). In the second algo- 
rithm a sequence of nonlinear operations was partly dis- 
entangling the state of the quantum computer by trans- 
formations of the type 
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|i> + |i>|i> 



(1) 



Nonlinear evolutions that lead (with an arbitrary accu- 
racy) to the required modification of entangled states can 
be obtained in a Weinberg-type nonlinear quantum me- 
chanics There is a problem, however. Using exactly 
the same trick [i.e. transformation it was shown in 
|l5| that (|l|) is responsible for arbitrarily fast influences 
between noninteracting systems (this should not be con- 
fused with another effect discussed by Gisin which 
was a result of the projection postulate). Therefore, the 
algorithm of Abrams and Lloyd makes an explicit use of 
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an unphysical and arbitrarily fast process, so there is a 
danger that this is the reason why it is fast. 

Fortunately, as we shall see below, this is not the case. 
To prove it we shall concentrate on the second algorithm. 
We will consider a concrete example of a nonlinear dy- 
namics and will use a formalism that is known to elim- 
inate the unphysical influences [^-21 1. A detailed dis- 
cussion of both algorithms in this context can be found 
in M. 



II. SECOND ABRAMS-LLOYD ALGORITHM 



Step 1. We begin with the state 

-- |0i,...,0„)|0) 



(2) 



where the first n qubits correspond to the input and the 
last qubit represents the output. 

Consider the unitary transformation acting as follows 
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c/|i) = -i=(-|o) + |i; 



step 2. 
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The input constists now of a uniform superposition of all 
the numbers < n < 2" — 1. 

Step 3. 

\m) = ^i^[i]> (7) 

where F is some unitary transformation (an oracle) that 
transforms the input into an output; /(ii, . . . , i„) equals 
1 or 0. 

Step 4- We assume that f{x) — 1 for at most one x. 
Denote by s, s = or 1, the number of x's that satisfy 
f{x) = 1. The state (||) can be written as 

|^[2]) = -^ E |0l,»2,...,Zn)|/(0l,l2,...,»n)) 

+ -4^ E 111, i2,.--,in)|/(li,i2, ■■•,««)) (9) 

12 ■■■in 

Let us note that with very high probabihty the state is 
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E (|0l,Z2,...,Z„)|0) + |ll,i2,...,«„)|0)) (10) 



With much smaller probability it is either 
1 



1 
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(|0i,Z2,...,Z„)|l) + |ll,i2,...,*n)|0)) (11) 



or 
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,i„)|0) -t- |ll,«2, 



(12) 



and is never in the form 
1 
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(|0i,Z2,...,z„)|l)-l-|li,i2,...,«„)|l)) (13) 



since this would mean there are two different numbers 
satisfying f{x) = 1 which contradicts our assumption. 
The idea of the algorithm is to apply to the flag qubit 
a nonlinearity that leaves (|l^) unchanged but (jT^) and 
( p^ ) transforms into 



^ |0i,i2, . . . ,l„)|l) + |ll,i2, 
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(14) 



Although such a dynamics can be approximated by a 
nonlinear Schrodinger dynamics of a Weinberg type, it is 
easy to show that it is unphysical |15|. 

To see this assume that the flag qubit does not interact 
with the n input ones. To simplify the discussion take 
n = 1 and consider the transformation (|l[). We assume 
that the nonlinear evolution is applied locally to the flag 
system. (^ implies the following transformation of the 
reduced density matrix of the flrst qubit 

1 



2(|0i)(0i| 



|li>(li 



^(|0i) + |li))((0i| + (li 



and therefore a fully mixed state evolves into a pure one. 
It can be shown IQ that Weinberg's description of sep- 
arated systems leads to 2-particle Schrodinger equations 
that induce this kind of behavior at a distance ( "faster- 
than-light telegraph"). 

Before performing a more detailed analysis let us il- 
lustrate the crucial element of the algorithm on a simple 
example. Take n = 3 and /(llO) = 1. The oracle pro- 
duces 

8-1/2 [ 
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The nonlinearity now "looks" at the second and the third 
input slots and sees the above kets as the following pairs 

8"^/2[|000)|0) + |100)|0) + 
|001)|0) + |101>|0) + 
|010)|0) + |110>|1) + 
|011)|0) + |111)|0)] (16) 

Now it scans each of the rows and does not do anything 
when two flag O's occur, but when it "notices" one and 
one 1 it changes to 1. So after this step we get 

8^1/2 [|000)|0) + |100)|0) + 
|001)|0) + |101)|0) + 
|010)|1) + |110)|1) + 
|011)|0) + |111)|0)] (17) 

Now the nonlinearity looks at the first and the third slots 
and sees the kets as the following pairs 

8-1/2 [|000)|0) ^ |oio)|l) + 
|001)|0) + |011)|0) + 
|100)|0) + |110)|1) + 
|101)|0) + |111)|0)] (18) 

It again behaves as before and what we get after this step 
looks as follows 

8-i/2^|000)|l) + |010>|1> + 
|001)|0) + |011>|0> + 
|100)|1) + |110>|1> + 
|101)|0) + |111>|0>] (19) 

Finally our nonlinearity looks at the first and the second 
slots and the state regroups in the following way 

8"i/2[|000)|l) + |001>|0> + 
|010)|1) + |011>|0> + 
|100)|1) + |101>|0> + 
|110)|1) + |111>|0>] (20) 

Now each row contains one 1 and in the final move all flag 
O's are switched to I's and the state partly disentangles: 

8-1/2 [|000) ^ |QQ^^ ^ 

|010) + |011) + 

|100) + |101) + 

|110) + |111)]|1) (21) 

Of course, in case s = the entire state does not change 
during the operation and a measurement on the flag qubit 
gives with certainty. Such an algorithm is fast and 



allows to distinguish between s = and s ^ in a linear 
time. The number of operations is of the order of n as 
compared to 2" typical of a slow data-search algorithm. 

It follows that we have an algorithm that is fast but si- 
multaneously makes an explicit use of an arbitrarily fast, 
physically unacaptable process. Now I will show that the 
faster-than-light effect can be eliminated without any loss 
in the efficiency of the algorithm. 

To do so I will use an explicit nonlinear dynamics and 
apply it locally to the flag qubit. By saying that the dy- 
namics is applied locally it is meant that we are using 
an appropriate (n -f l)-particle extension of a nonlinear 
1-particle dynamics. Assume the n + 1 subsystems do 
not interact with one another. The extension is local 
if it satisfies the following condition: A reduced den- 
sity matrix of any of the n + 1 subsystems satisfies a 
Liouville-von Neumann (nonlinear) equation which con- 
tains the reduced density matrix and Hamiltonian of only 
this particular subsystem. In the case of Weinberg's non- 
linear quantum mechanics of pure states and for finite- 
dimensional systems the extension of this type was in- 
troduced by Polchinski [[l7| . Its generalization to mixed 
states was given in [ [L8[|19| |. The extension to more gen- 
eral theories was discussed in [2^j2^ and applied to con- 
crete problems in |^,^. The algebraic origin of the 
Polchinski-Bona- Jordan formulation was discussed in de- 
tail in and It should be noted that extensions 
of this type are nonuniquc if one starts with a nonlin- 
ear dynamics of state vectors (cf. |ljj2^ ) . They become 
unique if the dynamics is from the outset given in terms 
of density matrices. From the point of view of the al- 
gorithm the uniqueness problem is irrelevant so we can 
stick to the nonlinear Schrodinger equation framework. 



III. MODIFICATION OF THE FOURTH STEP 

The above procedure can be, in principle, implemented 
in terms of a very complicated Schrodinger-type dynam- 
ics of the entire {n + l)-particle system of the quantum 
computer. It cannot be achieved by applying the nonlin- 
ear evolution locally to the flag qubit without generating 
the unphysical influences between different parts of the 
computer. 

Below I propose a simpler procedure which is based on 
a nonlinearity which is applied only to the flag qubit. 

Let us begin with a 1-qubit system whose dynamics is 
described by the nonlinear Schrodinger equation 

i\^) = etanh(^Q!(?/;|A- r/l|V'))A|V') (22) 

where e is the magnitude of the nonlinearity, a is a very 
large real number (say, a ~ 2"), 

A = r^{\0){0\ - |1)(1|) + yT3^(|0)(l| + |1)(0|), (23) 

and T] is small but nonzero. For lip) = |0) the expression 
under tanh vanishes. For a small admixture of |1) and 
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sufficiently large a the mobility with a nonzero frequency 
begins and an arbitrarily small amount of |1) can be suf- 
ficiently amplified. The Polchinski-type local extension 
of the dynamics to the entire quantum computer is ||26|] 

i\i>) = etanh ® (A - yyl)]*) j ® 

The {n + l)-particle solution is 

|*f> = cos - (g) Asinwi^ |*o> (24) 



with 



= etanh ^aTrp(A — ijl) 
-^tanh(iS) 



(25) 



where 



p==Tri...„|*o)(*o| 



2" 



2" 



is the reduced density matrix of the flag system after the 
first three steps of the original Abrams-Lloyd algorithm. 
The average of (73 = |0)(0| — |l)(f | at the flag subsystem 



(a3) = (vl'4|l(")®a3|vl/,) 

2""^ — s 2"^"'" — s 
= — cos2wt + 2?7^ —ain^ ujt (26) 

For s — the average is constant in time and equals 1. 
For s ~ 1, rj^ ^ 0, and sufficiently large a it oscillates 
with w w e. For t Ri ir/e the average is (as) « —1, which 
means that almost all flag O's in (^) have been changed to 
I's. Therefore instead of applying a complicated nonlin- 
ear dynamics of the original "Step 4" one can use the fact 
that in the local description the 1-particle nonlinearity is 
sensitive to the reduced density matrix of the particle. 

This kind of algorithm cannot distinguish between dif- 
ferent nonzero values of s, but clearly distinguishes be- 
tween s = and s 7^ in a way that is insensitive to 
small fluctuations of the parameters. It is interesting 
that the modified "Step 4" is essentially non-algorithmic 
and more resembles an effect of puncturing a baloon than 
performing an algorithm. 
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